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The electromagnetic response of bilayer graphene in a magnetic field is studied in comparison with 
that of monolayer graphene. Both types of graphene turn out to be qualitatively quite similar in 
dielectric and screening characteristics, especially those deriving from vacuum fluctuations, but the 
effect is generally much more sizable for bilayers. The presence of the zero-(energy-)mode Landau 
levels is a feature specific to graphene. In bilayers, unlike in monolayers, the effect of the zero-mode 
levels becomes visible and even dominant in density response as an externally-controllable band gap 
develops. It is pointed out that the splitting of nearly-degenerate pseudo-zero-mode levels at each 
valley, specific to bilayer graphene, is controlled by an applied inplane electric field or by an injected 
current. In addition, a low-energy effective gauge theory of bilayer graphene is constructed. 
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I. INTRODUCTION 

Graphene, a monolayer of graphite, has recently been 
attracting great attention, both experimentally^^ and 
theoretically^i 5 -^^ for its unusual electronic transport, 
characteristic of "relativistic" charge carriers behaving 
like massless Dirac fermions. Graphene is naturally of 
interest from the viewpoint of relativistic quantum field 
theory, and is a special laboratory to test the particle- 
hole picture^ of the quantum vacuum and, especially in 
a magnetic field, to study peculiar quantum phenom- 
en a 10 ! 11 ' 1 ^ 13 ' 14 tied to the chiral and parity anomalies. 
Actually, the half-integer quantum Hall (QH) effect and 
the presence of the zero-energy Landau levels observed^ 
in graphene are a manifestation of spectral asymmetry 
implied by the anomaly. 

The bilayer (and multilayer) of graphene is as in- 
teresting and exoti o 15 i 16 as the monolayer. In bilayer 
graphene interlayer coupling modifies the intralayer rel- 
ativistic spectra to yield a quasiparticle spectrum with 
a parabolic energy dispersion^ The relativistic feature 
thereby disappears but the particle-hole structure still 
remains, leading to a "chiral" Schroedinger Hamilto- 
nian, which has no analog in particle physics. Bilayer 
graphene, like monolayers, is intrinsically a gapless semi- 
conductor but has a notable property that the energy 
gap between the conduction and valence bands is con- 
trollabl o 17 ' 18 ' 19 ' 20 ' 21 by use of external gates or chemical 
doping. 

Theoretical studies show that electronic transport 4 ^ 5 -! 6 - 
and screening 2 ^ in graphene are substantially different 
from those in standard planar systems. The difference 
becomes even more prominent under a magnetic fields 
For graphene the vacuum state is a dielectric medium 
with appreciable electric and magnetic susceptibilities 
over all range of wavelengths. Curiously the zero-energy 
Landau levels, though carrying normal Hall conductance 
e 2 /h per level, scarcely contribute to the dielectric effect. 

The purpose of this paper is to study the electromag- 



netic response of bilayer graphene in a magnetic field at 
integer filling factor v, in comparison with that of mono- 
layer graphene. It turns out that both types of graphene 
are qualitatively quite similar in dielectric and screening 
characteristics, but the effect is much more sizable for 
bilayers (because of the difference in the basic cyclotron 
energy). 

In bilayer graphene there arise two species of zero- 
energy levels (or zero modes), which, when an interlayer 
field is applied, move up or down (oppositely) at the 
two valleys. Remarkably, the dielectric effects due to 
the (pscudo) zero modes, unlike in monolayers, become 
visible and even grow steadily as the tunable band gap 
develops. It is pointed out that the splitting of nearly de- 
generate pseudo-zero-mode levels at each valley, specific 
to bilayer graphene, is controlled by an inplane electric 
field or by an injected current. In addition, we construct 
out of the response a low-energy effective gauge theory 
of bilayer graphene in a magnetic field and verify that 
the electric susceptibility of a QH system is generally ex- 
pressed as a ratio of the Hall conductance to the Landau 
gap. 

In Sec. II we briefly review the low-energy effective 
theory of bilayer graphene and study its Landau-level 
spectrum. In Sec. Ill we examine the electromagnetic re- 
sponse and screening properties of bilayer graphene. In 
Sec. IV we derive an effective gauge theory. In Sec. V 
we study the effect of an inplane field on the almost de- 
generate zero-mode levels at each valley and discuss its 
consequences. Section VI is devoted to a summary and 
discussion. 



II. BILAYER GRAPHENE 

The bilayer graphene consists of two coupled hexag- 
onal lattices of carbon atoms, arranged in Bernal A'B 
stacking, with inequivalent sites denoted as (A, B) in the 
bottom layer and (A',B ! ) in the top layer. The elec- 
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tron fields in the bilayer are described by four-component 
spinors on the four sites and, as in the case of monolayer 
graphene, their low-energy spectrum is governed by the 
electron states near the two inequivalent Fermi points K 
and K' in the Brillouin zone. The intralayer coupling 
70 = 7^4 B ~ 2.9 eV is related to the Fermi velocity in 
monolayer graphene, v = (V3/2) w 10 6 m/s, 

with the lattice constant cil = 0.246nm. The interlayer 
coupling 71 = 7^4/ b and 73 = ^ab 1 are one-order of mag- 
nitude weaker than the intralayer coupling 70; numeri- 
cally^ 7l w 0.30 eV and 73 « 0.10 eV. 

Actually the neighboring sites A 1 and B form interlayer 
dimers via 71 and get shifted to higher energy bands, 
and the low-energy sector is essentially described by two- 
component spinors residing on the A and B' sites. The 
effective Hamiltonian is written asi^ 
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where coupling to electromagnetic potentials (A t , A ) is 
introduced through H = H x — iH y , W = H x + iH y , and 
Hi = —idi + eAi. Here £ = 1 refers to the K valley with 
H+ = Kf=i and ^ = (V'AjV'b')* while £ = — 1 refers to 
the K' valley with H- = Wf=-i and x = (Xb'jXa)*- We 
suppress the electron spin, which is taken care of simply 
via spin degeneracy g s = 2. 

In Ho the first term with a linear dispersion represents 
direct interlayer hopping via 73, with characteristic veloc- 
ity ^3 = (v3/2) a^ 3 /h ~ wo/30, while the second term 
with a quadratic dispersion represents A B' hopping 
via the dimer state, giving h 2 /(2m*) — v 2 /"f\. 

The Has takes into account a possible asymmetry be- 
tween the two layers, leading to a gap U between the 
conduction and valence bands. An important feature of 
bilayer graphene is that such a gap is controllable 17 ' 18 ! 19 
by use of external gates, U rs eAAo with an interlayer 
voltage AA = A op - ^ ottoro . The second term in H &s 
is a layer asymmetry associated with the depleted charge 
on the AB' dimer sites; we call it a kinetic asymmetry. 

Let us place graphene in a strong uniform magnetic 
field B > normal to the sample plane; to this end we set 
Ai(x) — > A B = B (—y, 0). It is convenient to rescale II = 
VteB a = (y/2/i) a and II* = (V2/£) at with the mag- 
netic length t = \j\feB so that \a,a [ \ = 1. The kinetic 
terms thereby acquire the scale U3 — > v^^/2/£ = uj^ and 
l/(2m*) -> eB/TO* = w c ; numerically, w c w 3.9 x B[T] 
meV and W3 » 1.2 x ^/£?[T] meV, where B[T] stands for 
the magnetic field in Tesla. The associated Landau-level 
spectra scale like io^y/n and uj c yjn(n — 1) with the level 
index n = 0, 1, • • •, with ratio ~ 0.3/(\/n~ _ T v /B[T]). 
Therefore the L03 term is practically negligible, compared 



with the u> c term, for higher Landau levels n > 2 under 
a strong magnetic field. With this in mind we cast the 
Hamiltonian TL^ in the form 
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with A = £w 3 /w c « ±0.3/^B[T]. Here fi = £L//(2w c ) 
stands for half of the band gap in units of uj c , and it 
seems feasible^ to achieve an interlayer- voltage change of 
magnitude fi ~ 0(1). Note that the kinetic asymmtry ~ 
/izata is very weak, with z = 2uj c /ji w 0.026 x B[T] <C 
1. 

Let us for the moment set a tiny parameter z — > 0. It 
is, then, generally seen from the structure of TL^ that its 
eigenmodes are the same as those of H $ l^^o and that the 
spectrum of H.^ is symmetric about e = 0, except for a 
possible e = w c /i spectrum or the zero-energy spectrum 
of Ti^l^o- We shall call the e = ±w c |/x| eigenmodes of 
TL^ pseudo-zero-modes, or simply "zero" modes. 

Such spectra are slightly modulated by the tiny O(^iz) 
kinetic asymmetry. For A = it is possible to write down 
the eigenmodes explicitly. The "nonzero" modes of 
are Landau levels of energy 
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labeled by integers n = ±2, ±3,..., and p x (or yo = 

£ 2 Px); A™ = m{i - (M - Hcrc s n = sgn{«} 

■ 1 specifies the sign of the energy e n . The associated 
eigenmodes are written as 

1 / \* 

^W,( X ) = ~J= {Cn<P\n\{x),-S n C-(j> lnl _ 2 {x:)) , (2.4) 



where = y/l ± /t„/e^ and = s„ y^fr 



1) 



A^; 

n (x) = (f)n{y — yo) (e 11 */' /\ / 2tt£ 2 ) are the eigenfunc- 
tions for the Landau levels (n, yo) of the usual Hall elec- 
tron. 

In bilayer graphene there arise two nearly-degenerate 
zero-mode levels per valley and spin, with spectrum 



u c fi, ei = w c fi(l - z), 



(2.5) 



and eigenfunctions ipoy a = ( < / , o( x )>0) t and ipi Vo — 
(0i(x),O)*. We take, without loss of generality, U/uj c = 
2£/i > and label the Landau levels associated with 
the £ = 1 valley (i.e., the -0 sector with fj, > 0) by 
n = 0+, 1, ±2, • • • and those associated with the £ = — 1 
valley (the x sector with fi < 0) by n = 0_, — 1,±2, • • •. 
Note that the zero-mode spectrum is ordered according 
to eo_ < £-i < < ei < e 0+ ; see Fig. 1. An interlayer 
voltage U oc 2fi thus works to shift the zero modes op- 
positely at the two valleys, opening a gap U, while the 
nonzero-mode levels get shifted only slightly and remain 
symmetric, apart from the O(^iz) asymmetry. 

The zero-mode levels at each valley are degenerate for 
z — > and the degeneracy remains even in the presence 
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FIG. 1: (Color online) Landau-level spectrum. The zero- 
mode levels move almost linearly with the interlayer voltage 
oc n and oppositely at the two valleys while other levels are 
shifted only slightly. 



of the linear kinetic term. This is because the zero modes 
persist even for A ^ 0; indeed, D(j> = with D — \(r — a 2 
has two independent solutions of power series in A, 



00' (x) = (1/No) 2_j a « 03n(x), 
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with normalization factors Nq — J2 n a n an d N 2 — 
SnCn! (n is given by a„ with replacement 2/3 — > 4/3 
and 1/3 — > 2/3 in the argument of the gamma functions. 
These n = 0',, 1' solutions have an infinite radius of con- 
vergence in A. This verifies that the index of the Dirac 
Hamiltonian, Index [7i^| M ^o] = dim KerD — dim KerZ?t, 
is two, apart from the Landau-level degeneracy, 



Index[H 5 | M ^ ] = 2 x 
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The zero-mode spectrum, corrected by the O(fiz) 
asymmetry, eventually reads 



e < = oJ c p{^ - 260(A)}, 
e v = uj c p{l - z6i(A)}, 

with 60(A) = J2n 3na l/J2„ a l and hW 
l)^/£„C>r 

60(A) = (1/2)A 2 + 0.05A 4 + ---, 
61(A) = 1 + (1/2) A 2 - 0.036 A 4 + 
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There is a common level shift of 0(zA 2 ) while the tiny 
level splitting ~ uj c pz <C oj c is practically unchanged (for 
A ~ 0.3). It is thus generally difficult to resolve the 



almost degenerate zero-mode levels; in a sense, it is the 
nonzero index (|2.7|) of the Hamiltonian H.^ that underlies 
this stability in the zero-mode degeneracy. 

In connection with the index, it would be worth re- 
marking that D'cf) = has no solution; a solution of 
infinite power series in 1/A, starting with </>o( x ), fails to 
converge for finite A. Actually the A — > 00 limit cor- 
responds to the case of monolayer graphene with linear 
dispersion. The index thus jumps from "2" to "-1" as 
one passes from finite A to infinite A. 

The linear kinetic term affects the spectrum e n of the 
\n\ > 2 levels only to 0(A 2 ), which is still negligible in 
a strong magnetic field. We shall therefore set A — » in 
most of our analysis below, except in Sec. V where we 
discuss possible resolution of the zero-mode levels. 

For actual calculations it is useful to make the Landau- 
level structure explicit via the expansio n 25 ! 26 ^(x, t) = 
L™,t,„( x Kyo) ^n(yo,t)- (From now on, we shall only 
display the tp sector since the x sector is obtained by 
reversing the signs of p and A.) The Hamiltonian H 
thereby is rewritten as 



H = I dy ^2tple n ip n , 
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and the charge density p~ P (t) = J d 2 xe lp ' x ip'tp as 



fc,n— — 00 

with the coefficient matrix 
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r = (ri,r 2 ) = (i£ 2 d/dyo,yo) stands for the center coor- 
dinate with uncertainty [ri,r2] = if 2 - Here 
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for k > n, and f n k{p) = [A-«(-p)] t ; P = PyUPxl actually 
fkn (p) are the coefficient functions for the charge density 
of the ordinary Hall electrons. 

In a similar fashion one can derive the expression for 
the current operator j ~ SH/SA. We omit it here, and 
simply remark that the current has a component coming 
from the O(pz) asymmetry as well. 



III. ELECTROMAGNETIC RESPONSE 

In this section we study the electromagnetic response 
of bilayer graphene. Let us first consider the polarization 
function P(p, a;)(~ — i (pp)) in Fourier space, 
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where = £fe — the sum is taken over occupied 
levels {n} and unoccupied levels {k}. In what follows 
we focus on the real part of P(p, uj) in the static limit 
lu —>■ 0, and denote the components coming from the 
virtual (n — ► k — > n) transitions as 
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Actually |5fc„(p)| 2 = ffnfc(-p) 9kn(p) are functions of 
£ 2 p 2 and are thus symmetric in (k, n), which implies the 
relation P% — —P^- 

For conventional QH systems the polarization function 
vanishes for the vacuum since the charge operator triv- 
ially annihilates the vacuum p\v = 0) = 0. For graphene, 
in contrast, p\v = 0) ^ because of pair creation from 
the Dirac sea, and even the vacuum acquires nonzero po- 
larization 
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p( P , o)u =0 = EE p -(p)- 



(3.3) 



look into the leading long- wavelength part ~ 0(p 2 ) of 
P(p,0), to be denoted as P( 2 )(p,0), which is related to 
the electric susceptibility a c — — (e 2 /p 2 ) P( 2 )(p, 0). A 
look into the matrix elements in Eq. (|2.12[) shows that 

(2) 

P„ (p, 0) derives only from virtual transitions to the ad- 
jacent levels (n — ► n ± 1) and the related ones across the 
Dirac sea (n — » — n± 1). A direct calculation then yields 
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for 2 < Inl < iV — 1, and 
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(3.7) 

(3.8) 

The bottom of the Dirac sea yields P^ ] N = (P (2) )^ J + 
(p(2))-(^- x ) oc AT - 1 + 0(1/N), which properly makes 

p (2) i,= = E"; " 

electric susceptibility (per valley and spin) 



P! (2) +Pr = -(p72™ c )2|/*| 



n=2 Pi„ finite. This leads to the vacuum 



Some care is needed in carrying out sums over an infinite 
number of Landau levels, which are potentially singular. 
For regularization we here choose, as before,— to trun- 
cate the spectrum to a finite interval —N<n,k<N 
and let N — > oo in the very end. Regularization not 
only keeps calculations under control but also refines 
them: Actually, instead of Eq. (|3.3p . one may first con- 
sider, for a given level n, the polarization "per level" 
P„ = J2k=~N Pn by summing over all levels k but n, 
and then obtain P(p, 0) = J2 n Pn by summing over filled 
levels n. This gives the same result as Eq. (|3.3|) owing to 
the antisymmetry P% = — PJ}. 

Equation (|3.3p is an expression for the tp sector (P — ► 
P^). For the charge operator p x = x^X m the x sector 
one may replace, in Eq. (|2.11[) . gkn(p) by 



9kJP) = 9-k,-n(p), 
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i.e., with the sign of p reversed. One can define (P x )^(p) 
accordingly and write the vacuum polarization function 

as P*(P,0)Uo = Ef =2 Elo( pX )"«(P)- Note that 
these components enjoy the following property 
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This implies, in particular, that the tp and x sec- 
tors contribute equally to the vacuum polarization 
P^(P,0)|„ =0 = P x (p, 0)1^0, and also that P(p,0) = 
P^(p)+P x (p) is the same for the charge-conjugate states 
with v = ± integer. (In view of this, we shall focus on 
the case v > from now on.) 

In calculating the density response we suppose that the 
almost degenerate zero-mode levels at each valley remain 
practically inseparable and treat them as both occupied 
or empty; accordingly we set z — ► below. Let us first 
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An alternative expression for this a™ c is obtained by 
summing up the virtual (— n —►nil) processes, or 
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where e' n = y n(n — 1) + u 2 . 

Equations (13. 6p and (|3.8p tell us that the susceptibility 
carried by a positive-energy level of the same n is slightly 
different by 0{p) terms at the two valleys. Such O(p) 
corrections are visible only for the zero-mode levels since 
the two valleys are practically indistinguishable for n > 2 
(up to O(pz) splitting). Let us suppose that the electrons 
fill up an integral number v of Landau levels, with uni- 
form density (p) = p = v/(2tt(. 2 ). We write v — J2 n v n m 
terms of the filling factors v n of the nth level [ < v n < 4 
for n > 2 and < ^{1,0} ^ 4; ^{1.0} = ^i+^o+Ij with both 
valley and spin taken into account. The electric suscep- 
tibility at integer filling v is then written as 
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FIG. 2: (Color online) (a) Electric susceptibility a c decreases 
with the band gap 2|/i| cu c for v = and rises for v =4, 8 
and 12. (b) Static susceptibility function a c [p], in units of 
e 2 /(2nuj c ), at v = 0,4,8, 12 for zero band gap oc /i = (thick 
curves) and for a finite gap — 0.5 (thin red curves), (c) 
Static dielectric function e(p, 0) — 1, in units of \/2a./{e^ui c t), 
at v = ~ 12 for a band gap oc /i = and — 0.5. 



(This a e is even in is, and applies to the case of holes 
with v < equally well.) 

It is amusing to note here that some expressions greatly 
simplify for a special value fi — 1/2: The spectrum is 
equally spaced, e„ = s n o; c (|n| — 1/2) for \n\ > 2 and 
e\ = £o+ = w c /2, apart from the 0(/j,z) splitting. For the 
susceptibility one finds (3 n — 1 for n > 2, and F(l/2) — 
1/2. This yields 



*e|/i=i/2 = (25i/q + ^)e 2 /(27rw c ), 



(3.13) 



which rises linearly with i/ = 4, 8, 12, • • -. 

In Fig. 2 (a) we plot cv e as a function of the band gap 
2|/x| w c for v — 0, 4, 8 and 12. The vacuum susceptibility 
a™ c = a c |„=o is almost comparable to the contribution 



of the filled n = 2 level for \i ~ 0. The 



decreases 



gradually as the band gap oc [i develops; this is clear intu- 
itively. In contrast, when the zero-mode levels get filled, 
i.e., at v = 4, a c starts to grow almost linearly with 
and such characteristic behavior persists for higher v as 
well. This is somewhat unexpected but is easy to under- 
stand: The zero-mode levels move linearly with the band 



gap oc \fi\ while the n > 2 levels are shifted only slowly, 
as seen from Fig. 1. The growth of a c is therefore due 
to a decrease in activation gap from the zero-mode lev- 
els. This implies that the dielectric effect due to the zero 
modes, though negligible for <C 1, becomes dominant 
for large gaps \/j,\ ~ 1. 

This is also the case with the full expression for the 
polarization P(p, 0), which one can evaluate numerically 
using Eqs. ((2~T3>1) and In Fig. 2 (b) we plot the 

susceptibility function a c [p] = — (e 2 /p 2 ) P(p, 0) at u =0, 
4, 8, 12 for zero band gap and for a finite gap. There, 
as in the monolayer case, the zero-mode levels scarcely 
contribute to cv c [p] for zero gap, but one sees clearly that, 
as the band gap oc /i develops, they make a e [p] distinct 
between the vacuum and the v = 4 state. 

The susceptibility is also related to screening prop- 
erties of graphene. Let us now turn on the Coulomb 
interaction v = a/(eb|x|) or v p — 27ra/(eb|p|) with 
a = e 2 /(47reo) « 1/137 and the substrate dielectric con- 
stant £b, and study its effects in the random-phase ap- 
proximation (RPA) . The RPA dielectric function is writ- 
ten as2£ 



e(p,w) = 1 - v p P(p, 



(3.14) 



Figure 2 (c) shows the static function e(p, 0) — 1 for v = 0, 
4, 8, 12, plotted in units of \J~2ol / \e\,u) c i) . Note first that 
there is no screening at long distances, e(p, 0) — > 1 for 
p — > 0, as is typical of two-dimensional systems. As wave 
vector |p| is increased, e(p, 0) grows rapidly, becomes siz- 
able for |p|^ ~ 1, and then decreases only gradually for 
larger |p|. Such profiles of e(p, 0) and a e [p] in Fig. 2 are 
qualitatively quite similar to those of monolayer graphene 
studied earlier.— 

Still there are some clear differences: (1) Note that the 
basic Landau gap of bilayer graphene, ui c as 45 B[T] K, 
is about one-order of magnitude smaller than the mono- 
layer gap w™ ono w 400 y/B[T] K at B=1T, or 
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Numerically it turns out that, for both monolayer and bi- 
layer graphene, the vacuum susceptibility a™ c is around 
3 in units of e 2 / (27r uj c ) and the peak value of e(p, 0)1^=0— 
1 is around 1.8 in units of \^2a / (ebu> c £) . This actually 
means that e(p, 0) and a c [p] are numerically more siz- 
able for bilayers than for monolayers. In particular, for 
the bilayer the peak value of e(p, 0) in Fig. 2 (c) would 
range from e(p, 0) w 9.5 ~ 18 for ^ = 0^12, with the 
choice e b ~ 4 and B = IT, or V2a/(e h ujJ) » 5.1. (The 
value of the substrate dielectric constant £b depends on 
the structure of the sample; with SiC>2 on both sides of 
the bilayer, e.g., one can set eb ~ £si0 2 ~ 4.) This implies 
that the Coulomb interaction is very efficiently screened 
(and weakened) in bilayer graphene. 

(2) For graphene, unlike standard QH systems, even 
the v = vacuum state has an appreciable amount of 
polarization e(p, 0) — 1 over a wide range of wavelengths, 
which reflects the quantum fluctuations or "echoes" of 
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FIG. 3: (Color online) Exciton spectra at u—0 and 4 for \p\ = 
0.1, in units of a/thi- The spectra are reduced from the bare 
ones (dashed curves) via screening. At v — 4 there are two 
spectra with excitation gaps ~ u c ( y2 + p? =p |/x|). 



the Dirac sea in response to an applied field. The echoes 
are in a sense "harder" for monolayer graphene^ 3 - for 
which e(p, 0) appears almost constant over the wave- 
vector range in Fig. 2 (c). The rise of the peak val- 
ues of e(p, 0) — 1 with filling factor v is more prominent 
for monolayers than bilayers; see Eq. (|4.6[) . These fea- 
tures reflect the difference in the underlying Landau-level 
structures. 

(3) In bilayer graphene, unlike in monolayers, the ef- 
fects of the zero-mode Landau levels become visible as 
the band gap oc 2\p\ develops and this makes the v = 
and v = 4 states distinguishable in e(p, 0) and a e [p\. 

In the RPA the response function is written as 
-Prpa(p,^) = P(p, w )/e(p, <*>), f rom which one can de- 
rive^ the inter-Landau-level excitation spectra corrected 
by the Coulomb interaction. Isolating from P(p, ui) 
one of its poles at to ~ — e„ and setting e(p, u>) — 

1 — tip P(p. ui) — > fixes the pole position of Prpa(p, w), 
rpa 

k,n 



£fc - e« + Ae fc ,„(p) with 



Ae fe ,„(p) 



*|P| 



'fg(p), 



e h £ 2e(p,0) 
^g(p) =^2gnk(-p)gkn(p)e~ x /x, (3.16) 

where x = £ 2 p 2 /2. 

The vacuum state at v = supports excitons associ- 
ated with the n = —2 — * (1,0+) or (—1,0-) — > 2 transi- 
tions with the excitation gap e cj c (y2 + /i 2 + jj) and 



! / g (p)=4(l- M /V2 + M 2 )%), 



(3.17) 



where h(x) = e~ x {1 + — 3)}. The v — 4 state 

supports excitons with energy e w o; c (\/2 + /i 2 ± /z) and 

^g(p) = 2 (1 T m/v^ + £i 2 ) M 21 ); respectively; see Fig. 1. 
Similarly, at f = 47i (n > 2) the n — * n + 1 transition 

with a gap Acjc"'' ~ e«+i — e n gives rise to excitons with 
j/ g (0) = (y/n — 1 + \/n + l) 2 . As a result of screening, 
Aefc in (p) oc l/e(p, 0), one would observe prominent re- 
duction in magnitude of the exciton spectra; see Fig. 3. 



IV. EFFECTIVE GAUGE THEORY 

In this section we study low-energy response of bilayer 
graphene and construct an effective gauge theory. Let 
us consider the Hall conductance which is read from a 
response of the form ^a}i(A x dtAy — A y dtA x )\ it is calcu- 
lated from the current-current correlation function and 
from Berry's phase as well.— A direct calculation similar 
to the one in the monolayer case^ 3 - shows that cth — > e 2 £ 2 



per electron so that oh 



V(27Tfi) 



- /h per filled 



level. The result is independent of A and z. 

Care is needed to determine cth carried by the vacuum 
state. As before, we truncate the spectrum and find that 
the bottom of the Dirac sea contributes ~ — (N — l)e 2 /h. 
This yields the vacuum Hall conductance 



(4.1) 



per valley and spin. This implies that nonzero current 
and charge are induce d 12 i 13 ' 14 in the vacua at the two val- 
leys, but they combine to vanish in the vacuum, leaving 
no observable effect. In general, the Hall conductance 
<th is cast in the form of a spectral asymmetry and the 
vacuum Hall conductance cr^ ac is related to half of the 
index (|277|) . 

The long- wavelength response of bilayer graphene with 
uniform electron density p is now summarized by the 
Lagrangian 

1 



L A = peA ~e'£'p-e^ p A^A p 
+ ^a c Ejj - ia m {A 12 f. 



(4.2) 



Here A detects the charge density p, a c is the electric 
susceptibility in Eq. (|3.12|) and a m oc e 2 £ 2 p/m* stands 
for the magnetic susceptibility probed by a local variation 
A12 — d x A y — d y A x about B. 

This response is essentially the same as that for mono- 
layer graphene,— except for the values of a e and a m - Ac- 
cordingly the effective theory also takes the same form, 
i.e., a theory of a vector field b^ = (60,61,^2), with the 
Lagrangian to O(0 2 ), 



A u e^ x d„b x + I b + -L & M e» vX d v b 



( 2 



+ - 



1 



(6 



1 



feO ) 



2£ 2 puj cS 2 
and the effective cyclotron frequency 
uj cS = e 2 £ 2 p/a c = ui c g(v), 



2£ 2 p 
5bi2 v 5b\ 



(4.3) 



4FV) + 

^{1,0} (Ml + Sn>2 V n A 



-, (4.4) 



where b^ = d f Jj ll - d v b^ £ 2 p = v/(2ir) and v = Y, p n- 
An advantage of bosonizatio n 29 i 30 is that it allows one 
to handle the Coulomb interaction exactly; it is in- 
cluded in Eq. (|4.3[) with shorthand notation 5bi2 v 5bi2 = 
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/ d 2 y £612(2) u(x-y)£6i 2 (y) and <5&i 2 = 612 - P- (We 
have omitted the a m term from L c g [6] since the Coulomb 
interaction overtakes it at long wavelengths.) 

This effective Lagrangian not only reproduces the orig- 
inal response (|4.2[) . but also shows that the Coulomb in- 
teraction v p — 27ra/(eb|p|) substantially modifies the dis- 
persion of the cyclotron mode at long wavelengths p —> 0, 

w(p) ^uj cS + ^(£ 2 pv p + ---)p 2 , (4.5) 

where • • • involves a m /(a e u) e g), if a m is recovered. This 
excitation spectrum is in good agreement with the RPA 
result (|3.16p at long wavelengths; the filling factor v = 
2n£ 2 p in Eq. (|4~5|) corresponds to f g (0) in Eq. ([3T6]) . 

For graphene the Landau levels are not equally spaced 
and the excitation gaps depend on the level index n 
or v. At v = An (n > 1), e.g., the minimum gap is 

Aw'"' = (^n(n + 1) + p? — \/n(n — 1) + p 2 ) u) c , and 

w e fi in Eq. (14. 4|) represents such an activation gap. 

In particular, at v = 4 where the zero-mode levels are 

filled, the effective gap cx 5(4) = l/[F(p) + //] derives 

solely from the vacuum fluctuations; it deviates from the 

true gap Awe = (\/2 + p 2 — /1) lj c by about 20% for 

p = and by 5% or less for p > 0.4. At v — 8, w c ff 
(2) 

and Awe agree within 1.1% for all p. The agreement 
is almost exact for higher v. It is somewhat surprising 
that, as in the monolayer case, an effective theory con- 
structed from the long-wavelength response alone gives 
an excellent description of the excitation spectrum. Here 
we verify again from Eq. (|4.4[) that the susceptibility a e 
of a QH system is generally given by a ratio of the Hall 
conductance to the Landau gap^2 

a c = u{e 2 /h)/io cS ^an/^\ (4.6) 

This implies that a c depends on v and p nontrivially 
while it grows like a c cx v^l 2 for monolayer graphene. 

V. ZERO-MODE LANDAU LEVELS 

In this section we take a close look into the proper- 
ties of the zero-mode levels at each valley, with a tiny 
splitting Lu c pz <C uj c . Special care is needed in consider- 
ing the response from such an almost degenerate sector, 
since even a weak field applied as a probe may affect the 
true eigenmodes. The first task therefore is to resolve 
the degeneracy of the (0' + , 1') sector by diagonalizing the 
external probe plus the O(pz) asymmetry, 

SH = - Jd 2 xeA oP + Jdyo^AH^, (5.1) 

where A7i as = uj c pz dia,g[— a^a, aa']. We suppose that 
A$ is slowly varying in space and retain only terms up 
to first order in d x , i.e., E = — 8 x Aq. In practice, it is 
convenient, without loss of generality, to take Aq(x, t) — > 
A (y,t), or E || E y . 



We take the linear kinetic term cx A into account (since 
it may potentially be important at such a low energy 
scale). Within the (0+, 1') sector 5H effectively turns 3 ^. 
into the matrix Hamiltonian (in yo space) 

-e Mvo , t )- Kb+ (X) + ( K g b ;$ ^f^), (5.2) 

with £ = elE y /V2, k = (1/2) u c /i z; c(A) = 
En=o V^+lanCn/iNom). Here &±(A) = h(X)±b (X) 
in terms of 60 and b\ defined in Eq. (|2.9[) . Numerically, 

c(A) = l-3x 10~ 5 A 4 + ---, 
b+(A) = 1 + A 2 + 0.086 A 4 + ---, 
MA) = 1 -0.014 A 4 + ---. (5.3) 

The 0(A 4 ) corrections are practically negligible for A ~ 
0.3. We therefore set c(A) « 6_ (A) w 1 and may keep 
the 0(X 2 ) correction to b + (X). 

The Hamiltonian (|5.2p leads to the level splitting 
±V k 2 + E 2 with the new eigenmodes 0" and 1" related 
to 0^ and 1' via the unitary transformation -0q» = 

-00' cos6* — ^ysuvO, with tan2# = e£ E y /(^/2 k). Here 
we see that an inplane electric field enhances the split- 
ting of the zero-mode levels, with a gap 

Ae = Lo c ^J{p zf + 2e 2 £ 2 E 2 /u; 2 . (5.4) 

This shows, unexpectedly, that the zero-mode level gap 
is controllable with an injected Hall current. 

For a rough estimate let us note the following: If we 
take p ~ 0.05, the level gap is u c pz ~ 5 x 10 _3 meV at 
B = IT, with pz - 1.3 x 10- 3 B[T]; an applied field of 
strength En ~ lV/cm leads to a level gap of roughly the 
same magnitude, as seen from 

V2eim/uj c - 10~ 3 x £[V/cm]/.B[T] 3/2 . (5.5) 

One can calculate the spectral weight ~ |(0"|/O p |l")| 2 
and derive the susceptibility function associated with the 
1" — > 0" transition, 

(ac[p])l " = 27rAe l+£ 2 / K 2 ' (5 ' 6) 

with x = £ 2 p 2 /2 and spin degeneracy g s — 2. (Here we 
have set A — > for clarity; there is no appreciable change 
for A ~ 0.3.) This a c [p], with its scale set by the tiny 
gap Ae -C lu c , essentially governs the dielectric property 
of the v = 2 state. 

Figure 4 shows a c [p] and the associated dielectric func- 
tion e(p, 0) — 1 at v = v\" = 2. It is seen that both 
of them decrease rapidly as E|| becomes strong in the 

sequence z) = 0, 1, 2, 5; this is because 

the 1" and 0" modes are chosen so as to diagonalizc 
/ d 2 x A Q p (for k — > 0). In Fig. 4 both a e [p] and e(p, 0)-l 
are plotted in units l/\pz\ ~ (40/|/x|)(l/B[T]) times as 
large as those in Fig. 2. Numerically, for p ~ 0.05 the 
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12 3 4 




12 3 



FIG. 4: (Color online) a c [p] and e(p,0) — 1 at v = 2, in 
units of e 2 /(27rw c |/x|z) and y/2a/ (eb£ u c \^\z) , respectively, for 
inplane electric fields of strength \/2el\E,\^\/ (uj c \^l\z) — 0,1,2 
and 5. An inset depicts an enhancement of level splitting by 
an inplane field. 

level gap ~ 60 mK at B = IT would become as large as 
~ 300 mK (~ 6.2 GHz h) by an inplane field of E\\ ~ 5 
V/cm. This would lead to a e [0] ~ 14 and (the peak value 
of e(p, 0) — 1) ~ 30 in common units, i.e., one-order of 
magnitude larger than those in Fig. 2. This suggests that 
the dielectric effect would show a marked enhancement 
around v ~ 2 if the level gap is resolved by an inplane 
field at very low temperatures. 

Actually, a larger gap, which leads to better level reso- 
lution, tends to suppress the effect. Alternatively it will 
be more practical to observe the field-induced level gap 
Ae(Eii) via the QH effect with an injected current; one 
would be able to resolve the v — ±2 Hall plateaus us- 
ing a suitably strong current (that may even suppress 
the dielectric effect). It will also be possible to detect 
the level gap directly by micro-wave absorption or via 
conductance modulations 32 associated with it. 



VI. SUMMARY AND DISCUSSION 

In this paper we have studied the electromagnetic re- 
sponse of bilayer graphene in a magnetic field at integer 
filling factor i>, with emphasis on clarifying the similari- 
ties and differences in quantum features between mono- 
layers and bilayers. Bilayer graphene has a unique feature 
that its band gap is externally controllable; this makes 
bilayers richer in electronic properties. 

The particle-hole picture of the vacuum state is one of 
the basic features specific to graphene, both monolayers 
and multilayers, and is not shared with conventional QH 



systems. In graphene even the vacuum state responds to 
an external field and acts as a dielectric medium, with 
the Coulomb interaction being efficiently screened over 
the scale of the magnetic length. Graphene bilayers and 
monolayers are qualitatively quite similar in the vacuum 
dielectric characteristics (apart from some differences 
that reflect the underlying Landau-level structures), as 
we have seen in Sec. Ill, but numerically the dielectric 
effect (for the vacuum and for v ^ as well) is generally 
much more sizable for bilayers because of t he diff erence 
in the basic Landau gap, wf /u)f ono w 0.1y / S[fJ. 

The presence of the zero-energy Landau levels is an- 
other basic feature specific to graphene. The monolayer 
supports four zero- mode levels (one per valley and spin), 
while the bilayer supports eight such levels (two per val- 
ley and spin). The zero- mode levels carry normal Hall 
conductance e 2 /h per level, but in monolayers their ef- 
fect is hardly visible in density response. In contrast, 
in bilayer graphene a gate-controlled interlayer field acts 
to open a band gap between the zero-mode levels at the 
two valleys and this valley gap, in a sense, activates them: 
The dielectric effect due to the pseudo-zero-mode levels 
grows linearly with the band gap and becomes dominant 
for large gaps, as we have seen in Sec. III. 

A finite band gap introduces an asymmetry in the zero- 
mode spectrum while it leaves other levels practically 
symmetric. The two zero-mode levels (per spin) at each 
valley thereby remain degenerate, apart from a tiny ki- 
netic splitting. This robustness in the zero-mode degen- 
eracy at each valley as well as the presence of the zero 
modes itself is a consequence of the nonzero index (j2.7[) of 
the basic bilayer Hamiltonian (with Aq — > and z — ► 0). 
In Sec. V we have pointed out that this tiny zero-mode 
level gap at each valley is enhanced (and even controlled) 
by an inplane electric field or by an injected current. Such 
a gap, if properly enhanced by an injected current, may 
be detected via the QH effect or directly via micro-wave 
absorption. 

Finally, in Sec. IV we have noted that the low-energy 
characteristics of bilayer graphene arc neatly summarized 
by an effective Chern-Simons gauge theory, which accom- 
modates graphene and standard QH systems equally well. 
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